The Effect of Interactions on the Conductance of Graphene Nanoribbons 
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We study the effects of the interaction between eiectrons and hoies on the conductance G of 
quasi-one-dimensional graphene systems. We first consider as a benchmark the limit in which 
all interactions are negligible, recovering the predictions of the tight-binding approximation for 
the spectrum of the system, and the well-known result G — 4 e 2 /h for the lowest conductance 
quantum. Then we consider an exactly solvable field theoretical model in which the electro-magnetic 
interactions are effectively local. Finally, we use the effective field theory formalism to develop an 
exactly solvable model in which we also include the effect of non-local interactions. We find that 
such interactions turn the nominally metallic armchair graphene nanoribbon into a semi-conductor, 
while the short-range interactions lead to a correction to the G = 4e 2 //i formula. 



I. INTRODUCTION 

Monolayer graphene is a truly two-dimensional (2D) 
system and a zero-gap semiconductor, where the elec- 
trons and holes behave as massless fermions with un- 
usual transport pr ope rties and display an anomalous 
quantum-Hall effect 1 2 . In order to utilize their impor- 
tant electrical features in certain applications, many ef- 
forts have been recently made to explore the proper- 
ties of lower-dimensional graphene nanostructures like 
nanoribbong2H2. By cutting graphene into a narrow rib- 
bon structure, the electrons and holes are laterally con- 
fined to form a quasi-one-dimensional (ID) structure, 
similar to the case of carbon nanotubes^and conventional 
semiconductor quantum wire d 8 " 13 '. The formation of ID 
subbands in graphene structures can lead to the quanti- 
zation of conductance. One of the most interesting ques- 
tions associated to this quantity is if the conductance G 
of the lower band deviates from that given by Landauer's 
formulas' for non-interacting systems, G — g s ^ where 
the degeneracy factor is g s — 4 for graphene and g s = 2 
for conventional systems. 

Effects of the electron-electron interaction in conven- 
tional ID systems were investigated by Tomonaga 1 -" 
many years ago, and since then a number of works have 
been published on this problemEa>22 One of the remark- 
able predictions of these theories is that, even in the ab- 
sence of scatterers such as impurities, the conductance 
G deviates from that given by Landauer's formula as 

2 

G = 75s X' wnere 7 depends on the interaction. In an 
experiment on a GaAs-AlGaAs quantum wire, Reilly et 
aZpShave found that 7 goes from 0.7 to 0.5 with increas- 
ing carrier density. Another intresting effect provided by 
Tomonaga's model is the possibility for electron-electron 
interaction to open a gap in the spectrum. This behavior 
has been recently experimentally verified in the contest 
of carbon nanotubes^S. 

Motivated by these experiments on quasi one- 
dimensional carbon-based systems, in this work we study 
the consequences of the electromagnetic interactions and 
quantum correlations between electrons and holes, inside 



nanoribbons built from a single layer of two-dimensional 
graphene. For these systems, tight binding calculations 
of the ID subbands lead to a dispersion relation 



E n (k) = ±hv^/k 2 + (n + a) 2 TT 2 /w 2 , 



(1) 



k is the longitudinal momentum, v ~ 10 6 to s _1 ~ 1/300 c 
is the band velocity and w is the nanoribbon width. 
n = 0, ±1,±2, .... is an integer for the subband index, 
and < \a\ < 0.5 depends on the crystallographic orien- 
tation of the GNR. The value of a — holds in armchair 
nanoribbons (AGNR) families with N=3p+2, where p is 
an integer and N is the number of dimer lines across the 
ribbon width. 

In the following, we mainly focus on the computation 
of the conductance of graphene quantum wires in the 
case in which n = and a ~ 0, so that, in the absence 
of interactions, the electron and hole n = bands form 
gapless Dirac cones. To this end, we develop a model 
motivated by an effective field theory (EFT) based on an 
expansion in kjki (or w/L where L is the length of the 
system), where fcy is the typical momentum of electrons 
and holes propagating in the transverse direction. 

We show that the combined effects of the interac- 
tion and of the low-dimensionality of the system drasti- 
cally change the structure of the ground-state and of the 
spectrum of excitations, which are typical of the two- 
dimensional graphene. In fact, the vacuum develops a 
condensate of electron-hole pairs and the spectrum of ex- 
citations turns out to be saturated by bosonic particle- 
hole collective modes. As a consequence, a gap opens 
up between valence and conduction band, turning the 
nanowire into a semiconductor: the AGNR behaves as 
a Mott insulator. Interactions at distances of the order 

of few lattice spacing lead to a correction to Landauer's 

2 
formula for the conductance, in the form G = 744- with 



T l + d/TT 



and — 7r < d < +00. 



The paper is organized as follows. In section [nj we in- 
troduce the model Hamiltonian for non-interacting elec- 
tron and holes and we calculate the free conductance 
of charge carriers. In section |IH| we develop a model 
in which the screened Coulomb interaction is effectively 



local, and compute the corresponding conductance. In 



section IV we use the effective field theory formalism to 
construct the most general Lagrangian describing the in- 
teraction at distance of electrons and holes in the wire. 
We show that, under suitable approximations, such a La- 
grangian reduces to an exactly solvable model, whose so- 
lution is given in section [V] The implications of these 
results on graphene physics will be presented in section 
VI Conclusions and perspective developments are sum- 
marized in section fVIII 



II. 



THEORY FOR GNR'S WITHOUT 
INTERACTION 



The role of dynamical electron-hole correlations in 
graphene quantum wires is expected to be very differ- 
ent than in mono-layer two-dimensional graphene. In 
fact, the free Dirac-like Hamiltonian for two-dimensional 
graphene has been shown to be a fixed point of the Renor- 
malization Group (RG) of the corresponding EFT^Zl. 
This implies that perturbative coupling constants of such 
an EFT can be made arbitrarily small by applying RG 
transformation. Hence, the free Dirac Hamiltonian can 
be used as a good starting approximation for discussing 
the dynamics of electrons and holes in a graphene layer, 
with interactions providing at most, small logarithmic 
corrections. On the other hand, irPSwe have shown that 
the lowest-order interaction terms the EFT describing 
quasi one-dimensional graphene have dimensionless cou- 
pling constants, and therefore appear at the same order 
in the EFT expansion of the free kinetic terms. This 
implies that interactions can never be neglected, and in 
general are expected to shape the physics of such quan- 
tum wires. 

In view of such considerations it is instructive to first 
consider as a bench mark the quantum motion near the 
K point of the electrons in the conductance band and 
holes in the valence band of GNR's in the limit in which 
all interactions are neglected. Such an analysis allows to 
clarify which properties of the wire are induced by the 
dynamical correlations. 

The second-quantized Hamiltonian leading to the de- 
sired free single-particle spectrum of the lowest energy 
band 



M 2 v 2 

u(k) = ±hv\/k 2 + — 2— 



(2) 



is of course 



Ho = vTi I dx ijj'(t,x) I - 



Mv\ 
io\d x + — I tp(t,x), (3) 



where ip(t, x) is the fermion field operator. In the follow- 
ing, we shall work in a natural system of units for this 
problem, in which K = v = 1. Note that, in such units, 
the speed of light is c = 1/(3 ~ 300. 

In order to exploit the formal analogy with the rela- 
tivistic Dirac theory, it is convenient to introduce position 



contravariant vectors 

x 11 = (vt,x) = (t,x), 
momentum contravariant vectors 

r = (E/v,p) = (E,p), 



(4) 



(5) 



and the metric tensor as g^ v — diag[l,— 1]. In addition, 
let 7 M and 7 s be 2 x 2 matrices obeying the usual Dirac 
algebra: 



{7^,7"} = 2gT 

{7 s , V} = 



(6) 

(7) 



and, (7 ) 2 = 1. For example, one may choose a repre- 
sentation in which 



7 
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7 S = 7°7 1 = en 



(8) 



e^lv 



Note that with this choice one has 7^7' s 

Using such a set of definitions, the action associated 
to the free "Dirac" Hamiltonian pi) can be cast in the 
familiar form 



S = J d 2 xtP U~P-M\ i/j. 



(9) 



Furthermore, in the presence of a classical external 
field A^ , the action of free electrons and holes inside 
the system is given by 



5' 



d 2 x ip [ilp- Ivl\ ip, 



where the covariant derivative is defined as 
D^B^ + ieA^. 



(10) 



(11) 



The electric conductivity of the GNR described by the 



action ( 10 ) is given by 
1 



K[a(w )9 )] = 



E ext (uj,q) 



(12) 



where A® (x) — ($(£), 0) is the potential of a weak 

external electric field, E e . 

Applying the linear response theory, one immediately 
finds 



d 2 x e^ £ (j\x)) Ar = 1 (m 10 m A% xt (q), (13) 

where TL^ 11 (q) is the vacuum polarization tensor, defined 
as: 

iW(q) = Jd 2 i e^ (n\T[r(x)f(0)]\Q). (14) 



In the non interacting model, this matrix element can be 
computed exactly and reads 



iIP"(£) 



liu 






i(q) 



where 



I(q) 



2 f 1 X(l - x) 

~ Q J M 2 -x(l-x)<? 2 



(15) 



(16) 



This result can be used to readily obtain the conductiv- 
ity of the system. After restoring the appropriate powers 
of v and h, we find our final results: 



If 



a(u,q) 



^5 



uj/v 



irk (uj/v + q)(uj/v — q) 



Hi) 



27i. 



5{u/v - q)I(u/v), 



and by taking the Fourier transforrrp3 



If 



<j(lu,x) 



/U)X\ 



.- cos I — I I(oj/v) 



(17) 



(18) 



When the mass term M is zero and the electron and hole 
bands are gapless Dirac cones, as happens for armchair 
GNR's with N=3p+2, in the limit w -> 0, we get for 
the DC conductance G (defined as the u) — > limit of 



If 



z 



) at zero chemical potential: 



a- T . ( M ) 

This result must then be multiplied for the factor g s = 4, 
to account for the spin and sub-lattice degeneracy: 



G^G = g s -. (20) 

h 

When there is a band gap, i.e. M ^ 0, we expect that 
the same result holds when M is small with respect to 
the Coulomb energy necessary to induce transport in the 
system. Otherwise, the DC conductance at zero chemical 
potential is zero. 

This discussion can be straightforwardly generalized to 
the case of finite electron chemical potential. In general, 
in the formalism of quantum field theory, the transition 
to finite density is achieved by adding a term to the zero- 
density action in the form 



S[i/),ip] — ¥ S[ip,tp] — ifi Id xtp(x)"fo'il}(x 



(21) 



Hence, in the free-theory case, the finite-density effect is 
simply that of shifting the po component which enters 
in the free fermion propagator by the electron's chemical 
potential fj,: po — > pa + /x. As a consequence, the first 
branch-cut singularity of I{q 2 ) in the complex plane is 
translated to the point qo = M — /x, on the real axis. This 




FIG. 1: If the typical momenta exchanged by fermions in 
the wire is much smaller than the screening length of the 
electro-static interaction, the interaction can be replaced by 
an effective local vector-vector vertex. 



means that, even for arbitrarily small Coulomb energies, 
the regime of finite conductance is reached as soon as the 
electron density p, becomes equal or larger than the gap 
M, as expected. 

In the next section we will take into account the effect 
of the interaction between electrons and holes on G in 
the case of M = 0. 



III. A MODEL WITH LOCAL 
CURRENT-CURRENT INTERACTIONS 

Let us now begin our study of the effects of the inter- 
actions between electrons and the holes. In this section, 
we consider a model which emphasizes the consequence 
of the screening of the Coulomb force inside the wire, and 
completely neglects the interaction at distance. 

If the electro-static interaction between electrons and 
holes is short-ranged, then the very low-energy electro- 
dynamics can be described by an effective vector-vector 
interaction — see Fig. [T] 



dm = x (*7m*) (*7 m *) 



(22) 



where the coupling constant d is inversely proportional 
to the inverse of the screening mass m s , i.e. 



e 2 /ml 



(23) 



Notice the analogy with the Fermi theory for weak 
decays, in which the vector-boson mediated weak- 
interaction is replaced by an effective local axial-vector 
coupling. 

In our 1+1 dimensional effective theory, the choice of 
interaction (22) defines the massless Thirring Model, i.e. 



Sth = d x ipi 



2 {inisf) ■ 



(24) 



The calculation of the conductance of the wire in such a 
model can be carried out analytically, even in the non- 
perturbative regime. 

As shown irJ22H23| it [ s possible to explicitly reformu- 
late the massless Thirring Model in terms of a massless 



bosonic free theory. In the present context, the bosons of 
such a theory can be physically interpreted as a electron- 
hole bound states. The fact that the spectrum of the 
wire in the Thirring Model is gapless implies that short- 
ranged vector correlations cannot turn the metallic arm- 
chair nanoribbon into a semiconductor. 

Let us now compute the conductance of the wire . 
Following^ we bosonize the system and map the ( 24 ) 
into 



°-\ 



d 2 £ 



(W + l&x) 2 



l d + l -)(^r 



(25) 



where 9 is the scalar field obtained through the bosoniza- 
tion of the fermionic fields, while x ancl 4 1 are the bo- 
son fields introduced using the auxiliary vector field 
technique^. 

The bosonization dictionary for the massless Thirring 
modePSES states that the vector current is represented 
as follows 



J a = *7m* 



1 

27 



-\IV 



d v t 



(26) 



Hence, it is straightforward to calculate the response of 
the system to an external field, in linear response theory. 
The current-current correlation function reads 



interaction is kept finite. To this end, we rely on the ef- 
fective field theory (EFT) formalism, which represents a 
powerful tool to describe the low-energy dynamics of ar- 
bitrary closed systems — for a pedagogic introduction to 
EFT's see^, for a more technical treatment see e.gP2 — . 
Any time a physical system is characterized by a large 
separation in its relevant energy-momentum scales, the 
low-energy observables are expected to be insensitive to 
the details of the physics which involves only the hard 
degrees of freedom, with energy and momenta above the 
gap. In this case, the EFT formalism can be used to 
systematically construct a model Lagrangian which de- 
scribes the dynamics of the relevant low-energy degrees 
of freedom. In general, the EFT is much simpler to solve 
than the corresponding (more) fundamental theory, as it 
is formulated in terms of fewer degrees of freedom. In ad- 
dition, the advantage of the EFT description is that its 
Lagrangian can be systematically constructed, starti ng 
from symmetry arguments. Indeed, Weinberg theorerrPZl 
implies that one should consider the most general La- 
grangian, compatible with the symmetry properties of 
the underlying fundamental theory and with the funda- 
mental principles of quantum field theory: 



£ = J2 c >( A )d t 



(29) 



{J^{p)Jv{-p)) = [g^u 



V^Vu 



p 2 J 7r(l + d/n) ' 



(27) 



Using the procedure (13)-(18l we can compute the 



quantum conductance for the AGNR. The result is 
G = 9s 



h(l + d/n) 



(28) 



Notice that this results differs from Landauer's formula 
by the factor 1 1/(1 + d/%) which is due to the presence 
of short-ranged correlations. 



IV. INCLUDING THE EFFECTS OF THE 
INTERACTIONS AT DISTANCE 

In the previous section, we have developed a phe- 
nomenological model in which the Coulomb interaction 
inside the wire was assumed to be effectively zero-ranged. 
Such a model becomes reliable only if the typical mo- 
menta of the electrons and holes propagating in the wire 
are much smaller than the inverse screening length. 

We now want to improve on such an approach, by 
defining an alternative model in which the range of the 



In this Eq., Oi are field dependent local operators, A 
is an ultra-violet cut-off, chosen in the gap between the 
high-energy and the low-energy modes and Cj(A) are the 
(running) coupling constants. The price to pay in the 
EFT approach is that the Lagrangian contains in princi- 
ple an infinite number of new unknown parameters Cj (A) , 
which implicitly embody the information about the ultra- 
violet physics, above the cut-off. Such coefficients have 
to be calculated from the underlying microscopic theory, 
or have to be fitted from experiment. 



Clearly, the Lagrangian (29) does not yet represent a 



physical theory, since it depends on an infinite number of 
parameters. On the other hand, the presence of a large 
separation in the energy- momentum scales assures that, 
in order to compute observables to any finite accuracy, 
one needs to specify only a finite number of such effec- 
tive parameters 2 . Typically, this corresponds to the cou- 
pling of the lowest-dimensional operators Oi . In the case 
of perturbative theories, one can show that higher di- 
mensional interaction operators lead in general to higher 
order contributions 3 in the expansion of the observables 
in power of p/A. Hence, any desired accuracy can be 



1 This result holds in the chirally conserved phase, i.e. d > — n. 
We expect that a modification of Landauer's formula will exists 
even in the chirally broken phase—, but this scenario will not be 
discuss here. 



2 The situation is completely analog to the multi-pole expansion of 
classical electrodynamics: the electro-magnetic field at distances 
much larger than the size of the source can be determined with 
arbitrary accuracy from a finite number of multipole coefficients. 

3 At most, perturbative quantum fluctuations can lead to correc- 
tions to the naive dimensions of the fields. 



reached by considering effective Lagrangians with only a 
finite number of effective interactions. 

Let us now apply the EFT formalism to our specific 
case, in which the natural cut-off scale is provided by the 
transverse momentum in the wire, A — kx- We are inter- 
ested in an effective field theory describing only the dy- 
namics of low-energy degrees of freedom inside the wire, 
i.e. the electrons and holes propagating along the lon- 
gitudinal direction. Hence, in our approach, we propose 
a purely 1 + 1 dimensional description of the electronic 
properties of the wire, using the EFT technology to build 
an effective Lagrangian which simultaneously takes into 
account for the short- and long-ranged correlations in- 
duced by the presence of the transversal confinement and 
of the interactions between electrons inside the wire. 

The only way to include interactions at distance, while 
respecting the symmetry, causality, unitarity and local 
charge conservation constraints is to introduce an effec- 
tive gauge-boson field a^ , which we shall call the pseudo- 
photon, owing to the formal analogy with QED. It is 
important to stress that this is an effective field, which 
does not represent physical photons propagating inside 
the wire. Its role is to mediate the long distance interac- 
tions arising from the interplay of electro-magnetic cou- 
pling and transverse confinement, in a way that charge 
conservation is enforced locally, in the wire. In addition, 
since we are interested in the response to an external 
electric field, we include into our effective Lagrangian the 
coupling with a [physical) external electro-magnetic field 
A^ x . On the other hand, we neglect the dynamics of 
photon radiation, from the electrons in wire, as this a 
relativistic effect and arguably very small, in the limit of 
low electron density. 

The pseudo-photon field a p propagates at a speed v' = 
f IP' <l/0, which does not in general correspond to that 
of the electrons and holes. This is equivalent to say that 
we can define a gauge in which the free classical Euler- 
Lagrangian Eq.s for a M are: 



a 



df 2 



dx 2 



0. 



(30) 



This can be achieved by considering a free action in the 
form 



So[5„] = / d 2 x [ -- A f^r 



1 ; 
4 



where 



and 



J i-iu @[iQ'i' (Jv&fii 



a — O0) 

a\ = p' a\. 



(31) 

(32) 

(33) 

(34) 



The action (31) leads to the free wave Eq. pOt upon 



imposing the modified Lorentz gauge condition, 



Clearly, since the pseudo-photon and the massless 
fermions propagate at different velocities, the resulting 
theory will not be invariant under any "Lorentz" -like 
symmetry. 

The starting point to construct an effective Lagrangian 
which conserves locally the chage of the fermions is to in- 
troduce a gauge invariant coupling of the pseudo-photons 
and fermions: 



D„ 



df+ + ig a^ + ieA 



ext 



(36) 



Let us now perform a naive dimensional analysis to 
identify the lowest order operators in our EFT. In I + I 
dimensions, the fermion field has mass dimension 1/2, 
while the a M and A^ field have mass dimension I. 
The lowest-dimensional gauge invariant and " Lorentz" - 
invariant operators are therefore 

i>ifol>, m(A)tpip, {tpO^j) 2 , (0 = l,i'y s ,%,%'Ys) 

UvF v , TpidJ^-y^id^F^Tp. (37) 

The first, second and third lines contain operators of di- 
mension 2, 3 and 4, respectively. 

Form this point ahead we will consider a special config- 
uration for the AGNR, which geometrically assures the 
disappearance of the mass term, as in planar graphene, 
due to the contact between valence and conduction 
bands. As shown irPSthis hypotheses is fulfilled in 3p+2 
family of AGNR, that represents 1/3 of the whole possi- 
ble armchair nanoribbons. 

The role of the contact terms (VJO^) 2 is to mimic 
the ultra-violet physics which sets in when electrons and 
holes interact at a distance of the order of the inverse cut- 
off, i.e. of the transverse size of the wire, ~ 1/&t- The 
interaction between electrons and holes separated by dis- 
tances much larger than the inverse cut-off is mediated 
by the coupling with the pseudo-photon field, e.g. the 
ip Upij) term. 

The exact non-perturbative solution of the theory de- 
scribed by the EFT Lagrangian containing all the terms 
in (37) is of course a formidable problem. On the other 



hand, our main purpose is to study the combined effect 
of long-range and short-range correlations on the elec- 
tronic properties of the wire. From such a stand point, 
we choose to retain in our model Lagrangian only the 
vector- vector contact interaction term ('tpj^ip) 2 , the min- 
imal coupling term tp ilftp and the pseudo-photon kinetic 
term f^f^. 

In conclusion, the effective action of our model for the 
internal quantum electrodynamics of the wire reads 



d 2 i 



ip ilptp 



2 0H^) 2 



Z f ft™ 



(38) 



d^ = 0. 



(35) 



We emphasize the fact that such an action was obtained 
by retaining only some of the effective couplings. This 
choice unavoidably introduces some model dependence in 



our calculations. On the other hand, the gain is that the 
model defined by ( 38 ) is exactly and analytically solvable, 



even in the non-perturbative regime. In the following, we 
shall refer to the model define by the action (38) as to 
the Schwinger-Thirring model. 



V. EXACT NON-PERTURBATIVE SOLUTION 
OF THE SCHWINGER-THIRRING MODEL 

In this section, we discuss the exact analytic solution 4 
of the Schwinger-Thirring model — see also the discussion 
irPSL. 



A. Solution of the Schwinger Model 

Let us begin by reviewing the solution of the pure 
Schwinger model, i.e. for d = 0, which can be found in 
standard quantum field theory textbooks, such aip21. A 
remarkable feature of such a theory is that its classical ac- 
tion is invariant under two independent transformations, 
defined by the following different types of local rotations 
of the fermion fields 



i>(x) -> e ix{x) i>{x) ip{x) -> ^(x)e~ tx{x) 
%l>{x) -> e^^VOc) ${x) -> i>{x)e i ' lS ^ x 

and by the corresponding gauge transformations: 

a M -> cif, d^x, 

9 



e^PQ. 



(39) 
(40) 



(41) 
(42) 



The dynamical consequences of such a symmetry be- 
come evident once one parametrizes the photon field de- 
grees of freedom as 



9 V p 



(43) 



and re-expresses the path integral in terms of the fermion 
fields ip, V» and of the \ an d 4> fields. 

The invariance of the classical action defines a classi- 
cal symmetry of the system. Such a symmetry is realized 
also at the quantum level only if the full path integral 
remains unchanged under the transformations (39)- (40), 
and (41 1-(42 ). To verify if such a condition is realized, let 



us analyze the transformation properties of the fermionic 
measure T>ipT>ip. It is possible to show that the gauge 



symmetry ( 39 1 leaves invariant the fermionic measure 



In order to keep the notation as simple as possible, all the deriva- 
tions reported in this section correspond to the choice /3' = 1. 
The generalization to arbitrary values of /3' is straightforward, 
although the notation is considerably more involved. 



and therefore the symmetry defined by the transforma- 



tion (39l-(41) is respected also at the quantum level. On 



the other hand, under the chiral gauge transformation 
fl40|)-||42l one has 



VipVip -> VipVipJ- 



where 



J- 



exp 



1 



dfx -= &*4> mfd u 
g 2 



(44) 



(45) 



is the functional Jacobian determinant of the chiral trans- 



formation (40) and m = gl\f^ is the so-called Schwinger 
mass. Hence, the chiral gauge symmetry is said to be 
" anomalous" , i.e. broken at the quantum level. 

As a result of the chiral anomaly, the path integral of 
the Schwinger model can be written as 

Z = J D^Vi>V<f> exp [iS shw ty, #, <f\] , (46) 

where the Schwinger action <S' s j lw is defined as 



-Sshwl^'^0] 



d 2 i 



ipiqti]) 



1 1 

2 s 2 



= ~d»<P [d 2 -m 2 ] 9 M 



(47) 



As shown iiP2l, this partition function can be mapped 
into a purely bosonic one, which reads 



^shwl 6 *' 1 



d 2 i 



o^e 



,2 fl 2 



(48) 



Notice that this action is quadratic in the boson fields, 
so every correlation function can be computed exactly. 



B. Solution of the Schwinger-Thirring Model 



We are now in a condition to discuss the solution of 



our Schwinger-Thirring model ( 38 ) . Applying the chiral 
and scalar gauge transformations (39)-(41) and (40)- (42), 



and evaluating the corresponding anomalous shift in the 
fermionic measure, the action transforms into 



Sshw^'V^] 



Si 



ipi^ij) - 



- i W {d 2 



2 g 2 



2 (V'T/xV') - 
m 2 ) Bu <(> 



(49) 



As in the case of the Schwinger model, the strategy to 
solve this theory is to define a bosonization scheme such 
that the resulting formulation of the path integral is 
purely Gaussian. Such a procedure involves determin- 
ing the bosonized representation of the vector current, 
- 07 M - 0- The procedure is illustrated in the appendix (|Aj). 



The result of such an analysis is that under bosonization 
the vector current operator becomes 



*7m* 



1 



CavdyO- 



(50) 



i 



Having established the "bosonization dictionary'**^ 
for the fermionic current, it is immediate to obtain the 
bosonized version of the Schwinger-Thirring model. In- 
deed, Eq. (49) transforms into 



s s tie,cf>} = lJd 2 xUi + ^)(d fl e) 



, (51) 



with m 2 = g 2 /it. Formally, this action has the same 
form of the original Schwinger model bosonic action. The 
only difference is in the factor (1 + d/Tr), which rescales 
the kinetic term. This factor contributes to correlations 
functions and provides corrections to the mass of the 
Schwinger boson, but does not modify the qualitative 
structure of the ground-state and of the excitation spec- 
trum. 



PROPERTIES OF THE GNR IN THE 
SCHWINGER-THIRRING MODEL 



VI. 



In section |IV| we introduced the Schwinger-Thirring 
model for the electrodynamics of the graphene wire. In 
the previous section, we have shown that such a model 
can be mapped into a modified Schwinger model. The 
effect of the short-range vector-vector interaction which 
are present in our model is absorbed into the coeffi- 
cient (1 + d/ir), which multiplies the kinetic energy term. 
Hence, using the fact that the Schwinger model is ex- 
actly solvable, we are finally in a condition to predict 
some important properties for the wire. 



A. Ground-state structure and bosonization of the 
spectrum of excitations 

As in the original Schwinger model, the vacuum is char- 
acterized by a finite fermion condensate 



(ft|*tf|fi) t^O. 



(52) 



In the context of graphene nanowire physics, this re- 
sult implies that, even in the presence of short-distance 
vector-vector interaction, the electron-hole pair density 
in the vacuum is not zero. 

This effect is due to the anomalous breaking of the chi- 
ral symmetry, which is also responsible for the bosoniza- 
tion of the spectrum. The Schwinger-Thirring model 
contains no fermion excitations, but only an arbitrary 
number of free fermions-antifermions bound states, with 
mass M — g 2 /ir(\ + d/ir). In the context of graphene 
physics, this means that the spectrum of excitations of 



the nanowire starts with a single collective electrons-holes 
mode with dispersion relation 



w(k) = ±hv 




2 a/fi 



e^l + d/n) 



(53) 



and contains a continuos of multi-boson excitations, 
starting at the two-boson threshold. Additional thresh- 
olds for multi-boson excitations are located at nm, with 
n = 3, 4, 5... In these formulas, we have restored the con- 
stants v, c, h and a = ^j— — 1/137 is the fine structure 
constant. 



B. Conductance of the graphene wire 

The calculation of the conductance of the wire in linear 
response theory in the presence of electron-hole interac- 
tions can be performed following the same steps taken in 
the free case. It is important to recall that the model we 
are considering contains two distinct vector fields: the 
effective pseudo-photon field a M and the physical exter- 
nal photon field A^. The conductance is stimulated by 
the latter field, while the long-range dynamical correla- 
tions are generated by the former field. Consequently, 
the coupling of such two fields with the fermion field ip 
is parametrized by different coupling constants, g and e. 
Note that this is not the case in the original Schwinger 
model, in which one considers only one type of coupling 



of the fermions to the vector field — see e.g 



42 



As 



a consequence, the current-current correlation function 
which enters in the definition of the conductance in our 
model is not the same correlation function which enters 
in the Dyson series associated to the pseudo-photon mass 
renormalization 5 . 

The current-current correlation function IT™ of the 
Schwinger system to an external vectorial perturbation 
was calculated by Schwinger in his original worlP^ k u t, 
the bosonisation technique offers a different way to ob- 
tain the same result. By definition we have 



Z sh [A^ = 



i8Ap(x) ib~A v (y) 



A u =0 



logZ[A M ], (54) 



Va^VWD 1 ® e 



iS sch +i e ( d 2 x^{x)'y^(x)A^(x) 



(55) 



applying the bosonization technique, we find 

%n$(x,y) = --(e^eii^dreiy)). (56) 

n 



5 Note that, in our previous worH^SIthis point was overlooked. The 
calculation of conductance reported in such a work but holds only 
for values of the chemical potential larger than the Schwinger 
mass (i.e. of the dynamically generated semi-conductor's gap). 



After analytically continuing back to real time, the 
Fourier transform of such a result is 



„-rrSh/~N _ e " /-2, - - ' 



iH-nAq) = (Q 9iiv ~ QnQv 



7T q 2 — m 2 



(57) 



where m is the Schwinger mass. The effect of long-range 
correlations on current-current correlations can be read- 
off by comparing this formula its free counterpart, Eq. 



(151. 



By following exactly the same procedure one can com- 
pute the current-current correlation function in the com- 
plete Schwinger-Thirring model, i.e. using the action 
(pit. We find 



in?. (9) 



1 



7T f + d/lT 



{tg^u - q^u) 



1 



q 2 — rri 2 



(58) 



We note that the existence of a single pole implies that 
the conductance is entirely saturated by the propaga- 
tion of confined electron-hole bound states (Schwinger 
bosons), and not by (quasi) free electrons and holes. 

Since the system develops a gap, the conductance be- 
comes different from zero for electron densities larger 
than the gap, i.e. for /i > m. If is this the case, the 
procedure used from Eq. (15) to Eq. (19) can be re- 
peated to find the quantum of conductance. The result 



is 



G = g s 



h(l + d/n)' 



(59) 



where we have included the degeneracy factor. 

The present analysis has illustrated how the transport 
properties of the wire are determined by two different 
mechanisms which are related to the long- and short- 
range part of the interaction, respectively. Interactions at 
distance induce a gap in the spectrum, implying the insu- 
lating property of the nano-wire. However, long-ranged 
interactions alone are not able to modify the value of 
conductance's quantum. Corrections to the free-theory 
value 4e 2 //i value are completely due to short-ranged in- 
teractions. 



We have then taken into account the effect of the in- 
teraction by developing a model inspired by an effective 
field theory based on an expansion in k/kr, where k and 
kx are the momenta in the longitudinal and transverse 
directions, respectively. We have shown that long-range 
interactions dynamically generate a gap in the spectrum, 
turning the GNR into a Mott insulator. In addition, 
short-range interactions lead to a renormalisation of the 
free theory result for the conductance. Once both types 
of interactions are taken into account, one obtains that 
when the Fermi energy exceeds the dynamically gener- 
ated gap, the conductance is G — 4e 2 //i(l + d/tr), where 
d is the short-range interaction strength. 

A possible development of the present work would be 
to investigate how the properties of the wire change as 
a function of the temperature. Also, our field-theoretic 
approach can be easily implemented in the case of Carbon 
Nanotubes. 



Appendix A: Bosonization of the 
Schwinger-Thirring Lagrangian 

Let us consider the theory defined by the partition 
function 



^*,*,o„ I-^-mI 



P*P*X>a M exp 



d 2 x \ru v 



(Al) 



Notice that, for sake of simplicity and without loss of 
generality, we have dropped the "tilde" notation and we 
have assumed the same coupling constant e for both the 
dynamical a M and the external A^ fields. 

To solve this theory, we can proceed in the same way 
as section (V A), applying the gauge transformation also 



to the external field A^. Hence, we define: 
1 



A„ 



(d^a-ie^dvP). 



(A2) 
(A3) 



VII. CONCLUSIONS 

In this work, we have studied the effect of the interac- 
tions between electron and holes on the conductance G 
of quasi-one-dimensional graphene systems at zero tem- 
perature. We first considered the case in which all the 
interactions are absent and there is a gap in the electron 
and hole bands, described by an effective mass term M 
in the free Hamiltonian. In this case, when the Coulomb 
energy of the field applied to the ends of the wire is 
larger than the gap M, we recover the well known result 
G = 4e 2 /h. On the other hand, when such a Coulomb 
energy is much smaller than the gap, we obtain a finite 
conductance G = 4e 2 //i only when the chemical potential 
/i becomes equal or larger than M, as expected. 



The corresponding path integral becomes 



exp 



{ 



d 2 x^d^d 2 d^ 



+ !<*(- d^(x + a)- ie^du {4> + /?) )] *} 



(A4) 



In the presence of the external field, the application of the 
gauge and chiral gauge rotations leads to two anomalous 
terms, which are related to the <f> and j3 fields, respec- 
tively: 



h^ 2 



Z7T 



(A5) 



Consequently the path integral is written as 



1 



VdV<\> exp 
1 



d 2 x-(d,ef 



Z7T Ze Z7T 

Now we perform a shift in the held 0: 

8^6 + <J)/tt + /3/vr 



(A6) 
(A7) 

(A8) 



and we integrate over the held 0, as in the case of the 
pure Schwinger model. We obtain 



Z s 



= / VQt 



-Jd 2 ii(a M 8) 2 -^ + lm¥+iraV 2 



(A9) 



where m 2 = e 2 /ir. The term proportional to /3 2 can be 
dropped, as it does not contribute to correlation func- 
tions. The resulting partition function is 



Z e [ 



V6 e 



-/d ! ii(8,,fl) 2 + ira 2 8 2 - 



2 i i„2a2 a^fta^s 



(A10) 



We note that, in such a representation, the interac- 
tion with the external held is described by the term 
-j^d^fid^O. After re-expressing such a coupling in terms 
of the original external held A^ we find 



ZelA^} = / V9 



- J d 2 a;i(a M 9) 2 + im 2 e 2 -i-i=A„€„ fi a M e 



(All) 



This result shows that, after bosonization, the vector cur- 
rent operator becomes 



V-y^ 



j,vd v 6. 



(A12) 
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